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We present a complete phenomenological model accounting for the evolution of the cosmic-ray 
spectrum and composition with energy, based on the available data over the entire spectrum. We 
show that there is no need to postulate any additional component, other than one single Galactic 
component depending on rigidity alone, and one extragalactic component, whose characteristics 
are similar to those derived from a study of particle acceleration at mildly relativistic shocks in 
a GRB environment (Globus et al., 2015). In particular, we show that the resulting cosmic ray 
spectrum and composition satisfy the various constraints derived from the current data in the 
Galactic/extragalactic transition region, notably from the measurements of KASCADE Grande and 
Auger. Finally, we derive some generic features that a working phenomenological scenario may 
exhibit to give a global account of the cosmic ray data with a minimum number of free parameters. 


I. INTRODUCTION 

One century after the discovery of cosmic rays (CRs), 
their study remains one of the main focus of high en¬ 
ergy astrophysics and astroparticle physics. In the re¬ 
cent years, an important set of new cosmic-ray data 
have become accessible, thanks to major experimental 
progresses. In particular, the Pierre Auger Observatory 
(Auger) [I] and Telescope Array [2] experiments have ex¬ 
plored the ultra-high-energy cosmic rays (UHECRs) with 
unprecedented observational power. Below the ankle, 
KASCADE-Grande (KG) [3] has made crucial measure¬ 
ments of the cosmic-ray spectrum, including a separa¬ 
tion between low-mass and high-mass nuclei, which can 
be linked to measurements in the knee region made in 
particular by the KASCADEexperiment. Below the 
knee, new data have been available as well, notably from 
('ll HAM [a and TH ACEItjGj and most recently from 
PAMELA [Tj and AMS [5], allowing a more precise deter¬ 
mination of the relative abundances of the various nuclei 
among the Galactic cosmic rays (GCRs). 

These new data allow a more complete descrip¬ 
tion of the CR phenomenology over the entire energy 
range and in particular at the Galactic/extragalactic 
(GCR/EGCR) transition, taking into consideration both 
the CR energy spectrum and the composition. Some re¬ 
cent studies have claimed that the data could not be ac¬ 
counted for without invoking, in addition to the main 
GCR and extragalactic UHECRs, one or more additional 
(Galactic or extragalactic) components dominating the 
flux at intermediate energies, between the knee and the 
ankle (see for instance IMG]), and whose spectrum, com¬ 
position and flux normalisation could be adjusted at will 
to reconcile the observations with one’s preconception 
regarding the low-energy and/or high-energy CRs. In 
this Communication, we argue against such a necessity, 
and explicitly show that the global CR spectrum can be 


fully described in a natural way within a two-component 
model, one Galactic and one extragalactic, without the 
need of introducing a new, hypothetical component. 

We base our model on a very simple and generic de¬ 
scription of the GCR component, whose properties de¬ 
pend only on the rigidity of the particles (that is, for 
ultra-relativistic particles, the energy-to-charge ratio), 
and on an EGCR component that is directly borrowed 
from a previous work on particle acceleration at mildly 
relativistic internal shocks of GRBsjT3], which consis¬ 
tently predicts an evolution of the composition compati¬ 
ble with Auger data. We show that this model provides 
a fair account of all the significant features of the cosmic 
ray spectrum and composition from below the knee to 
the highest energies. 


II. MODEL 

A. The extragalactic component 

One of the most important results of Auger is the 
observation of a transition from a light-dominated to a 
heavy-dominated composition between a few 10 18 eV and 
a few 10 19 eV [HI 13], As we first proposed in m, 
this can be interpreted as a consequence of a low-energy 
cutoff of the protons in the sources, due to an intrin¬ 
sic limitation of the acceleration process. As long as 
the acceleration of particles is governed only by elec¬ 
tromagnetic processes, all nuclei behave in exactly the 
same way if they have the same magnetic rigidity. Nu¬ 
clei of charge Z are thus expected to reach an energy 
Emax(2) = Z x E max (p), where K max (p) is the maxi¬ 
mum proton energy. It is thus natural to expect that 
nuclei of higher and higher mass dominate the UHECR 
source composition above the maximum energy of the 
lower-mass nuclei. In more realistic situations, however, 
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FIG. 1. Left: EGCR flux as a function of energy for H, He and different ranges of nuclei (indicated by the charge labels), as 
predicted by our acceleration model in GRBs, adjusted to the Auger data. The figure was taken from m- Right: Effect on 
the propagated EGCR spectrum of different assumptions for the cosmological evolution of the source density, in (1 + z) a , with 
a = 2.1, 2.6, 3.0 or 3.5, as indicated. 


such a simple model should be amended to take into ac¬ 
count the energy losses and photo-dissociation processes 
that may occur in the acceleration site. In a recent study, 
we developed a numerical model for the acceleration of 
UHECRs in the mildly relativistic internal shocks of a 
GRB [13] . We showed that the relatively high density of 
energetic photons in the acceleration site leads to signif¬ 
icant photo-dissociation, which has two important con¬ 
sequences in the present context: i) the resulting maxi¬ 
mum energy of the nuclei is not strictly proportional to 
Z , but also reflects their photo-dissociation rate; and ii) 
the spectrum of the UHECRs eventually injected by the 
source into the intergalactic medium is close to a hard 
power-law (roughly in E~ x below E max ), but while all 
composed nuclei have essentially the same spectral in¬ 
dex, protons have a significantly steeper spectrum. This 
is due to the secondary neutrons, which are mostly pro¬ 
duced by photodisintegration processes during the accel¬ 
eration. Indeed, the charged particles mostly escape from 
the acceleration region in the weak scattering regime, i.e. 
at the highest energies. On the contrary, the secondary 
neutrons are not confined by the local magnetic fields, 
and thus escape with their production spectrum (which 
is similar to that of the nuclei at the shock), flowing freely 
out of the source before decaying into protons (see p[3j 
for more details). 

The model consistently predicts the shape of the spec¬ 
tra of individual nuclei, including their high-energy cutoff 
at the source. We then convoluted individual source in¬ 
jection over the GRB luminosity function and used our 
UHECR propagation code, taking into account energy 
losses, photo-dissociation and magnetic deflections[TT]. 
to derive the propagated spectrum which can be observed 
on Earth. The result is shown in Fig. [T] which is taken 
from da, with a comparison to the Auger data. The 
shaded area corresponds to the so-called cosmic variance 
and represent the expected range for the flux of the differ¬ 


ent nuclei, including 90% of independent realisations of 
the model. As can be seen, the propagated proton spec¬ 
trum is indeed much softer than that of the other nuclei. 
This model reproduces fairly well the overall spectrum, 
and shows a clear transition from a proton-dominated 
composition at the ankle to a Fe-dominated composition 
at the highest energies. 

In order to limit the number of free parameters, we use 
this model of the EGCR component in the present study. 
However, we allow for different assumptions regarding 
the cosmological evolution of the sources. Figure [l] shows 
the total spectrum obtained when the source density in¬ 
creases as function of redshift as (1 + z) a , with a = 2.1, 
2.6, 3.0 or 3.5. As expected, changing the cosmologi¬ 
cal evolution does not affect the high-energy part of the 
spectrum, since the contributing sources are all located at 
low redshifts, due to the GZK horizon effect. However, a 
stronger source evolution implies a larger contribution of 
the EGCR sources at low energy. Since the correspond¬ 
ing flux is dominated by protons, larger values of a result 
in larger contributions of extragalactic protons, which in¬ 
fluences the composition at the GCR/EGCR transition. 
From our calculations, we found that a relatively large 
evolution, with a between 3.0 and 3.5, provides the most 
striking agreement with the composition measurements 
over the whole energy range. In the following, we use 
a = 3.5, which is fully compatible with the observational 
constraints. In particular we verified that the gamma- 
ray emission resulting from the intergalactic showers as¬ 
sociated with the propagation of the EGCRs does not 
violate the measurements made by the Fermi telescope. 
We should however note that our results do not neces¬ 
sarily imply that the EGCR sources must have a strong 
cosmological evolution. Similar changes of the EGCR 
spectrum could also be obtained for instance by modify¬ 
ing the assumed source luminosity function or the turbu¬ 
lence structure at the shock. Moreover, the assumption 
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that GRBs are the sources of UHECRs is not critical for 
the success of the model. The key feature on which we 
rely here is the prediction of a softer spectrum for the 
protons, which can be expected in other cosmic accelera¬ 
tors as soon as they involve a significant amount of nuclei 
and the matter or radiation density is large enough in the 
source environment. 


B. The Galactic component 

For the GCR component, we assume that all nu¬ 
clei have the same rigidity spectrum, consisting of a 
broken power law with spectral index x below an en¬ 
ergy E hreak (Z) = Z x .E break (p), and a spectral index 
x + Ax above that energy, up to an exponential cutoff in 
exp( E f Z E max ) above E max . 

The slope x and the relative abundance of the vari¬ 
ous nuclei are simply adjusted on the most recent avail¬ 
able data at low energy[TB . Note that, we only use data 
above 300 GV, given the evidence found in PAMELA 
data(H] °f a change of slope of the spectrum of H and 
He below this rigidity (also confirmed by the most re¬ 
cent AMS data[2U]). Above 300 GV, we found that 
an index x = 2.67 provides a good fit of the data 
for all nuclei. In particular, TRAGER [2TU52] finds an 
index of x = 2.67 ± 0.08 for the combined elements 
heavier than He, while CREAM [24] find an index of 
2.66 ± 0.02 protons, and 2.66 ± 0.04 for elements heav¬ 
ier than He|25], There is admittedly some tension with 
the measurements of CREAM in the specific case of He 
nuclei [241. namely :r(He) = 2.58 ± 0.02. However, the 
most recent AMS02|2nj results are in good agreement 
with x(He) = 2.67. Therefore, instead of leaving the 
spectral index free for each nucleus, we decided to stick to 
the most natural assumption that all nuclei have the same 
spectrum in rigidity, and simply determine the relative 
abundances from the observational data. A combined 
fit of CREAM and AMS02 data is used to normalise the 
proton flux while only AMS02 is used for He nuclei, and a 
combined fit of CREAM, TRACER and ATIC-2[2>] data 
(whenever available) is used for heavier nuclei (CNO, Ne, 
Mg, Si, S, Ar, Ca and Fe). While some small corrections 
might be due to spallative processes (see e.g. m) , we 
note that the model is anyway robust with respect to the 
assumption of a single index x for all elements, since a 
change of 0.03 in the value of x for an element would 
result in a mere 23% difference of its relative abundance 
after three decades in energy, which we verified has very 
minor impact on our results. 

The change of slope Aa; (the knee) is a major feature 
in the GCR spectrum; its origin may be on the side of 
the acceleration process, for instance a consequence of a 
reduction of the number of sources contributing at higher 
and higher energy (see [33] and references therein). Or 
it may be a feature of the propagation of the GCRs in 
the interstellar medium, through a change of the diffu¬ 
sion regime, or the effect of a Galactic wind (e.g. 25 and 


references therein; see also [T2j for a recent account). It 
could also reflect some inhomogeneity in the GCR flux, 
associated with the granularity and/or intermittency of 
the sources. In the current modelling, we do not at¬ 
tempt to give any interpretation of this change of slope. 
Ax is adjusted in order to obtain a good fit of the spec¬ 
trum of individual nuclei (or classes of nuclei), as given by 

KASCADE[271 m and KASCADE-Grande[25H22]• We 
find values of Ax between ~0.3 and 0.5 coupled to values 
of Eb re ak(p) between ~2x 10 15 and 4 x 10 15 eV . 

Finally, the only remaining free parameter of the GCR 
component model is the energy scale, E max , of the high- 
energy cutoff, which we adjust to ensure a good fit of 
the data in the GCR/EGCR region. In particular values 
of E max between ~ 5 x 10 16 and 1.5 x 10 17 eV allow 
the heavy Galactic component to become very low to 
negligible at the ankle (which thus marks roughly the end 
of the GCR/EGCR transition |351I36| ) and the combined 
light Galactic and extragalactic component to produce an 
ankle around ~ 10 17 eV, as observed in KG data j30}l3l]. 


III. RESULTS 

In this section, we show the results obtained with the 
above model, after setting the values of the parameters 
to E hleak (p) = 10 15 - 5 eV, Ax = 0.45 and E max = 6 x 
10 16 eV, and combining the resulting Galactic component 
to the EGCR model using a = 3.5. 

The resulting spectrum is shown in Fig. [2] The 
left panel shows a global view of the CR energy spec¬ 
trum above 10 115 eV. In particular the assumed spec¬ 
tra of the 4 most abundant Galactic species (H, He, O 
and Fe) are compared with satellite and balloon borne 
measurements from various experiments f CREAM [24. 
[25], PAMELA [19], PAMELA-CALO[37], ATIC-2[26] and 
TRACER [22J [22]), while the total spectrum is compared 
below the knee with various balloon borne or ground 
based experiments (RUNJOB[38j, JACEE[39J, ATIC- 
1[40], ATlC-2[26] and Tibet-IIIjH]. Taking into account 
the existence of some discrepancies between the various 
experimental results, our Galactic component is in good 
agreement with data for both the total spectrum and the 
different components below the knee energy. The right 
panel shows a close up view at higher energy including 
the knee, the ankle and the highest energy regions. Our 
model predictions are compared with KASCADE[28] as 
reconstructed with the QGSJetII-3 model [42] , KG[2I] as 
reconstructed with QGSJetII-4]221 and Auger [44] data. 
The dashed line shows the flux obtained by summing the 
H and He spectra, including both GCRs and EGCRs. 
This is in good agreement with the KG data, when 
only the so-called “light” component (corresponding to 
electron-poor showers) is taken into account, which is 
intended to correspond to H and He nuclei according to 
the analysis presented in G2B m- The dotted-dashed line 
represents the flux of all the elements heavier than Mg, 
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FIG. 2. Comparison between the spectrum of our global CR model and the available data, detailed in the upper right 
corners. The individual spectra of H, He, O and Fe are shown for the GCR component (left, with the indicated rescaling). In 
addition, the “heavy” and “light” components determined by KASCADE-Grande are shown with dashed and dotted-dashed 
line (see text). Right: close up view, showing the sum of the H and He fluxes (dashed) and the elements heavier than Mg 
(dotted-dashed). 


and is to be compared with the so-called “heavy” com¬ 
ponent. As can be seen, both components are consistent 
with the data, and fit well within the KG systematic er¬ 
rors. In particular, we find that the heavy component 
exhibits a knee slightly below 10 17 eV, while the light 
component shows an ankle slightly above that energy as 
the light EGCR component become dominant over the 
sharply decreasing light Galactic component. The be¬ 
haviors of both the heavy and light components are in 
very good agreement with the KG findings [251132) . 

Of course, the quantitative agreement between our 
model and the data depends on the underlying assump¬ 
tion regarding the hadronic model used to reconstruct 
KG data. However, the choice of the hadronic model 
mostly influences the relative normalization of the dif¬ 
ferent components. In particular, for all the hadronic 
models tested in pH [32], the heavy knee and light ankle 
features remain, and the post knee shapes of the light and 
heavy component remain similar to one another. More¬ 
over, we found that the relative abundances for these 
components predicted by our model agree best with the 
data reconstructed with the QGSJetII-4 which is, to date, 
the most recent of the models tested against KG data 
and the only one that includes recent experimental con¬ 
straints for LHC dat;0 

In addition to these spectral features, our model pro¬ 


1 EPOS-LHC, the most recent version of the EPOSHA. l46l model 
(including constraints from LHC data) has not yet been con¬ 
fronted to KG data. Given the quoted similarity between EPOS- 
LHC and QGSJetII-4 in terms of the predicted correlation be¬ 
tween the shower size and the muon number in KG energy range, 
one can however anticipate that the light and heavy components 


vides a detailed description of the CR composition from 
the knee to the highest energies, which allows a compar¬ 
ison with the data. The left panel of Fig. [3] shows the 
relative abundance of H, He and the following two domi¬ 
nant sets of nuclei, namely CNO (and sub-CNO) nuclei, 
and sub-Fe (and Fe) nuclei. For protons and sub-Fe nu¬ 
clei, we also show separately the Galactic component, 
using dotted lines. 

As can be seen, even though the Galactic protons es¬ 
sentially disappear at ~ 1.5 10 1 ' eV, the abundance of 
protons never drops below 15%, and rises up again to 
more than 50% (with a maximum around 60%) just above 
5 10 17 eV. The fact that these protons, which ensure a 
dominantly light component across the ankle, are extra- 
galactic protons, is fully consistent with the anisotropy 
measurement of Auger. Indeed, a Galactic component 
of protons would most probably produce a significant 
anisotropy towards the Galactic center and/or plane, 
which is excluded by the data. Finally, the proton frac¬ 
tion is seen to rapidly decrease above the ankle, to finally 
vanish above a few 10 19 eV, letting heavier and heavier 
nuclei dominate the UHECR spectrum. 

The behaviour of Fe (and sub-Fe) nuclei is quite differ¬ 
ent, as there is practically no overlap between the Galac¬ 
tic component, which ends at a few 10 18 eV (i.e. 26 times 
higher in energy than the Galactic protons), and the 


deduced from KG data should be quite similar for both models. 
Likewise these models should yield a lighter composition at the 
knee than previously estimated from KASCADE data|271 1281 , 
probably more consistent with direct measurements at lower en¬ 
ergies. 
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FIG. 3. Left: Relative abundance of H, He and the elements in the indicated charge ranges, as a function of energy. For 
H and the elements between Sc and Fe, the purely Galactic component is shown with faint dotted lines. Center: comparison 
between the model predictions for the evolution with energy of the depth of the shower maximum, X max , and the Auger data, 
for three different hadronic models. Right: same as central panel, for the X max variance, <7(X max ). 


extragalactic component, which rises up strongly above 
10 19 eV, to reach 60% at 10 2 ° eV. 


It is interesting to note that, according to our model, 
the dominant class of nuclei over roughly one decade in 
energy, between ~ 6 10 18 eV and ~ 5 10 19 eV, should in 
fact be CNO. This appears in very good agreement with 
the recent Auger findings m- 


The spectra of individual nuclei are unfortunately very 
difficult to measure separately, which prevents a direct 
comparison with the data. However, it is possible to 
compare the data with the model predictions for the 
composition-dependent observables, namely the depth 
of the maximum shower development, traditionally re¬ 
ferred to as A' max , and its spread (among the whole set 
of showers) at a given energy, cr(A' max ). This is done in 
Fig-! where we plotted the evolution of these two ob¬ 
servables (central and right panels) with energy, together 
with the Auger data. For this, we simulated the devel¬ 
opment of a large number of cosmic-ray showers for the 
different nuclei and energies, using the CONEX shower 
simulator |47j with three different choices of the hadronic 
interaction model (SIBYLL2.1 FT8] . QGSJetII-4[33] and 
EPOS-LHC[J5] 0S]. The agreement between the predic¬ 
tion of our model and the data is remarkable over the en¬ 
tire energy range, both qualitatively and quantitatively, 
especially when the shower development is calculated us¬ 
ing the EPOS-LHC hadronic model. It is again interest¬ 
ing to note that this model takes into account the recent 
constraints from measurements performed at LHC. Al¬ 
though they probably do not reproduce perfectly all air 
showers properties [49], the most recent hadronic models 
seem to give a more coherent picture of the evolution 
of the composition deduced from indirect measurements, 
from the knee to the highest energies. 


IV. SUMMARY 

We showed that the whole CR spectrum, including 
the key region of the GCR/EGCR transition, can be de¬ 
scribed by simply superposing a rigidity dependent GCR 
component and a generic EGCR model, without addi¬ 
tional degrees of freedom. 

In our model, the GCR component is identical for all 
nuclei with the same rigidity. The maximum energy of 
protons accelerated in Galactic sources is ~ 6 10 16 eV, 
and the transition towards extragalactic protons takes 
place around 10 17 eV, where KASCADE-Grande ob¬ 
serves an ankle in the light CR component. While the 
knee-like break in the GCR proton component occurs at 
~ 3 10 15 eV, the corresponding break in the Fe compo¬ 
nents appears at ~ 8 10 16 eV, which is in agreement with 
the observed “heavy-knee” in the KASCADE-Grande 
data. The normalisations of the light and heavy com¬ 
ponents are also in good agreement with the data. 

Our results suggest that extragalactic protons account 
for more than 50% of the total flux from ~ 5 10 17 eV 
to ~ 5 10 18 eV, and drop below 10% above 3 10 19 eV. 
The dominant class of nuclei between ~ 6 10 18 eV and 
~ 5 10 19 eV is CNO. The evolution of the composition 
predicted by our model has been shown to be fully com¬ 
patible with the Auger data fCT ITS] . across the observed 
transition from a light-dominated to a heavy-dominated 
composition. 

An important reason for the success of the model is 
the fact that the EGCR source spectrum is significantly 
steeper for protons than for the heavier nuclei. As re¬ 
called above, this is because most of EGCR protons in¬ 
jected in the intergalactic medium below ~ 10 19 eV, 
are in fact decay products of freely escaping secondary 
neutrons, produced during the acceleration through the 
photo-dissociation of heavier nuclei. While this is a direct 
consequence of our particle acceleration model, presented 
in detail in El, we believe that it is a generic feature of 
UHECR acceleration processes occurring in photon-rich 
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environments. As a consequence, the agreement between 
the model predictions and the data should not be taken 
per se as an argument in favour of the GRB source model, 
but in favour of the global CR interpretation scheme de¬ 
veloped in this study. Combined with a relatively low 
maximum proton energy, a single EGCR component can 
thus explain in a consistent way the evolution of the com¬ 
position above 10 18 eV as well as the spectral and com¬ 
position features observed below the ankle, down to the 
knee. In particular, the relative energy scale between the 
“light ankle” and the ankle (a factor ~ 30 ~ Z(Fe)/Z(p)) 
is most naturally explained within this framework involv¬ 
ing a transition between two components. 

We thus conclude that, unless new data will contradict 
the current observational status, there is no need for in¬ 
voking any additional CR component. Now, if one opts 
indeed for the minimal assumption that a single com¬ 
ponent accounts for all GCRs, then our results strongly 
suggest that this component should be able to accelerate 


protons up to at least ~ 5 — 610 16 eV. This might be in 
tension with the generally accepted models for particle 
acceleration at the shocks of individual supernova rem¬ 
nants in the Galaxy, which hardly reach energies much 
larger than 10 15 eV. This could suggest that these mod¬ 
els should be modified, or that other classes of sources 
may be have a dominant contribution to the GCRs, per¬ 
haps through collective acceleration processes as might 
be expected in superbubbles (see e.g. m and references 
therein). 
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